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I Abstract. By means of an effective relativistic nuclear equation of state in the 

framework of the nonextensive statistical mechanics, characterized by power-law 
quantum distributions, we study the phase transition from hadronic matter to quark- 
gluon plasma at finite temperature and baryon density. The analysis is performed by 
requiring the Gibbs conditions on the global conservation of baryon number, electric 
charge fraction and zero net strangeness. We show that nonextensive statistical effects 
" strongly influence the strangeness production during the pure hadronic phase and the 

hadron-quark-gluon mixed phase transition, also for small deviations from the standard 
i ^ '_ Boltzmann-Gibbs statistics. 

> 

C\| ■ 1. Introduction 

m 

One of the major challenges in high energy heavy ion collisions is the determination 
of the physical proprieties of strongly interacting nuclear matter at high temperature 
and baryon density. Relativistic heavy-ion collisions provide the unique possibility to 
explore in laboratory nuclear matter under extreme regimes in which the baryon density 
can reach values of a few times the saturation nuclear density and/or high temperatures. 
In these conditions, phase transition phenomena in the hot and dense fireball created 
during the collisions can take place [I]. Various results from QCD inspired models 
indicate that, increasing the baryon chemical potential in the phase diagram, a region 
of non-singular but rapid cross-over of thermodynamic observable around a quasi-critical 
temperature, leads to a critical endpoint (CEP), beyond which the system shows a first 
order phase transition from confined to deconfined matter. The existence or exclusion of 
a CEP has not yet been confirmed by QCD lattice simulations. Actually, there are some 
extrapolation techniques to finite chemical potentials, although the precise location of 
the CEP is still a matter of debate [2j [3]. Such a CEP can be in principle detected 
in future high-energy compressed nuclear matter experiments such as FAIR at GSI in 
Darmstadt p[] and NICA at JINR in Dubna [5j. In this direction interesting results 
have been obtained at low SPS energy and are foreseen at a low-energy scan at RHIC 

The process of deconfinement and the equation of state (EOS) of hot and dense 
nuclear matter can in principle be described by QCD, however such a theory is highly 
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non-perturbative in the energy density range involved in relativistic heavy-ion collisions. 
The generated quark-gluon plasma (QGP) in the early stages of the collisions does 
not at all resemble a quasi-ideal gas of quarks and gluons because strongly dynamical 
correlations are present, including long-range interactions [91 Q2JJ [11] . Therefore, in the 
absence of a converging method to approach QCD at finite density one has often to resort 
to effective and phenomenological model investigations to obtain qualitative results. 

In the last years there is a growing interest to high energy physics applications of 
nonextensive statistics pEJ E2 EH] and several authors have outlined the possibility that 
experimental observations in relativistic heavy ion collisions can reflect nonextensive 
statistical behaviors (T5J M, El CEE1 HSJ EQl EI] • 

Nonextensive statistical mechanics introduced by Tsallis consists in a generalization 
of the common Boltzmann-Gibbs statistical mechanics and it is based upon the 
introduction of entropy [H] 

S *W = i 1 ~ J [ /(X) ] ? dQ ) ' (/ /(X) dQ = X ) ' (1) 

where /(x) stands for a normalized probability distribution, x and dQ denoting, 
respectively, a generic point and the volume element in the corresponding phase space 
(here and in the following we set the Boltzmann and the Planck constant equal to 
unity). The real parameter q determines the degree of non-additivity exhibited by the 
entropy form ([Q) which reduces to the standard Boltzmann-Gibbs entropy in the limit 
q — > 1. By means of maximizing the entropy S q , under appropriate constraints, it is 
possible to obtain a probability distribution (or particle distribution) which generalized, 
in the classical limit, the Maxwell-Boltzmann distribution. The nonextensive classical 
distribution can be seen as a superposition of Boltzmann one with different temperatures 
which have a mean value corresponding to the temperature appearing in the Tsallis 
distribution [20^ 122] . 

The existence of nonextensive statistical effects incorporates in an effective way 
the presence of strong dynamical correlations, long-range interactions and microscopic 
memory effects and it affects significantly the finite temperature and density nuclear 
EOS (23] EH [25] [26] [27] EH] - In fact, by varying temperature and density, the EOS reflects 
in terms of the macroscopic thermodynamical variables the microscopic interactions of 
the different phases of nuclear matter. 

In a previous work, we have studied the relevance of nonextensive statistical effects 
in the hadron to quark-gluon phase transition by means of a relativistic nuclear equation 
of state in a range of finite temperature and density for which the production of 
strange particle can be negligible and we have studied the deconfinement transition from 
nucleonic matter to up and down quark matter [29]. By increasing the temperature of 
the system the producing of strangeness through associated production becomes much 
more important and we expect that nonextensive statistical effect considerably influence 
the strange particle densities in the EOS [30J. 

The main goal of this paper is to study how nonextensive statistical effects influence, 
from a phenomenological point of view, the nuclear EOS in presence of strange matter 
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and, as a consequence, the relative hadron to quark-gluon phase transition at finite 
temperature and baryon density reachable in high-energy heavy-ion collisions. At this 
scope we consider an effective relativistic mean-field model with the inclusion of the 
full octet of baryons, the A-isobar degrees of freedom and the lightest pseudoscalar 
and vector mesons. These last particles are considered in the so-called one-body 
contribution, taking into account their effective chemical potentials depending on the 
self-consistent interaction between baryons [3T| 132]. 

During the phase transition, the Gibbs conditions for the phase equilibrium are 
applied by requiring the global conservation of three "charges": baryon number, electric 
charge (isospin) and strangeness number. A crucial feature of such conservation laws 
is that, during the mixed phase transition, the baryon, electric charge and strangeness 
densities can be different in the single hadron and quark phase, although the total net 
charge densities are kept fixed [33] ■ As a consequence, the pressure in the mixed phase 
is in general not constant and the nuclear incompressibility does not vanish at all |34j . 

The nonextensive index q is considered here as a free parameter, even if, in principle, 
it should depends on the physical condition generated in the reaction, on the fluctuation 
of the temperature and it should be related to microscopic quantities such as, for 
example, the mean interparticle interaction length, the screening length and the collision 
frequency into the parton plasma. Moreover, let us observe that when the entropic 
q parameter is smaller than one, the high energy tail of the particle distribution is 
depleted; when q is greater than one, the energy tail is enhanced. Hence the nonextensive 
statistics entails a sensible difference of the power-law particle distribution shape in the 
high energy region with respect to the standard statistics. In the present work we 
will focus our study for small deviations from the standard statistics and for values of 
q > 1, because these values were obtained in several phenomeno logical studies in high 
energy heavy ion collisions (see, for example, Ref.s [201 I2U [35]). In this context, let us 
remember that, in the diffusional approximation, a value q > 1 implies the presence of 
a super diffusion among the constituent particles (the mean square displacement obeys 
to a power law behavior (x 2 ) oc t a , with a > 1) [36]. 

The paper is organized as follows. In Sections |2] and El we present, respectively, the 
main features of the effective nonextensive hadronic and quark-gluon EOS. In Section HI 
we report the Gibbs conditions applied for the hadron to quark-gluon phase transition 
and the consequent formation of a mixed phase. Finally, in Sections |5] and |6l we present 
the main results and summarize our conclusions. 



2. Nonextensive hadronic equation of state 

The relativistic mean-field model is widely successful used for describing the properties 
of finite nuclei as well as hot and dense nuclear matter [371 EH Ell HQ] . In this framework, 
the Lagrangian density describing hadronic matter can be written as 

£ = ^octct + £a + £qpm j (2) 
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where £ oc tet stands for the full octet of the lightest baryons (p, n, A, £ + , E°, £~, H°, E~~) 
interacting with cr, u, p meson fields; £a corresponds to the degrees of freedom for the 
A-isobars (A ++ , A + , A , A~) and £ qpm is related to a gas of the lightest pseudoscalar 
and vector mesons with an effective chemical potentials (see below for details). 
The Lagrangian for the self-interacting octet of baryons can be written as jlO] 

£octct = 7m ^ - ( M k - 9ak o-) - g uk 7m ^ ~ 9 P k 7m t' P M ] 4'k 

k 

+ g {d^a - m l° 2 ) - 3 o (9*n °? ~\ h (d.N cr 4 ) 
+ 2 m l + 4 c (91n u^) 2 + - m 2 p p M • p " 

- \F, U F^ - \G, V & V , (3) 

where the sum runs over the full octet of baryons, is the vacuum baryon mass 
of index k, the quantity t denotes the isospin operator that acts on the baryon and 
the field strength tensors for the vector mesons are given by the usual expressions 
Fp.v = d^ujy — d v uj^ G^ u = d^pu — d u pn. 

In regime of finite values of temperature and density, a state of high density 
resonance matter may be formed and the A(1232)-isobar degrees of freedom are expected 
to play a central role [4T1 1421 1431 144] . In particular, the formation of resonances 
matter contributes essentially to baryon stopping, hadronic flow effects and enhanced 
strangeness [45] . The Lagrangian density concerning the A-isobars can be expressed as 

HUE] 

£ A = ^a v - ( m a - 9a^<y) - q^ai^Wa , ( 4 ) 

where is the Rarita-Schwinger spinor for the A-baryon. Due to the uncertainty on 
the meson-A coupling constants, we limit ourselves to consider only the coupling with 
the cr and u> meson fields (g a A and 5Wi), more of which are explored in the literature 

[171 SHI SH]. 

The field equations in a mean field approximation are 

(i 7M 3" - M* - 9ua u> - g P ii%iP)4> k = , (5) 
(z 7 m^ -Ml- g uA j°u))ipA = , (6) 
m 2 a a + agl N a 2 + bg* N a 3 = ^ g ai p\ , (7) 

i 

m lu + cgtiN^ 3 = 9^ipf , (8) 



m 2 p 



9/nT3ipf , (9) 



where cr = (cr), a; = (u>°), p = (p°) and 5 = (63), are the nonvanishing expectation 
values of meson fields and the effective mass of the z-th baryon is defined as 

M* = Mi - g ai a . (10) 
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The pf and pf are the baryon density and the baryon scalar density, respectively. They 
are given by [50] 



/d 3 k 
——[n^-fiiik)}, (11) 

/c\ 3 k M* 
(2^8^-[nK*)+W)l- (12) 

where ji = 2Jj + 1 is the degeneracy spin factor (■joctet = 2 for the baryon octet 
and 7a = 4) and rii(k) and rii(k) are the g-deformed fermion particle and antiparticle 
distributions. Following Ref. [26], we can write, for f3(E*(k) — p*) > 0, 

ni{k)= [l + (q-l)P(ET(k)-rt)]W*-V + l' (13) 
and, for P(E*{k) -//*)< 0, 

nt(k) = [l + {l-q)ftE;(k)-ri)]W-ti + l ■ (M) 
The corresponding antiparticle distribution rij can be obtained with the substitution 
p* — > —p*- The baryon effective energy is defined as E*{k) = \fk? + M* A and the 
effective particle chemical potentials p* are given in terms of the meson fields as follows 

A =lk- 9u,iU - Qpihip , (15) 

where pi are the thermodynamical chemical potentials pi = de/dpi and % is the third 
component of the isospin of i-th baryon. 

Because of we are going to describe a finite temperature and density nuclear matter 
with respect to strong interaction, we have to require the conservation of three "charges" : 
baryon number (B), electric charge (C) and strangeness number (S). For this reason, 
the system is described by three independent chemical potentials: ps, pc and ps, 
respectively, the baryon, the electric charge and the strangeness chemical potentials. 
Therefore, the chemical potential of particle of index i can be written as 

Pi = bip B + Cipc + Sips • (16) 

where b i: Ci and Sj are, respectively, the baryon, the electric charge and the strangeness 
quantum numbers of the z-th hadronic species. 

The meson fields are obtained as a solution of the field equations in the mean field 
approximation and the related meson-nucleon couplings constant (g a N, g^N and g P N) 
will be fixed to the parameters set marked as GM3 of Ref. [10]. The implementation of 
the model with the inclusion of hyperons, require the determination of the corresponding 
meson- hyperon coupling constant, that have to be fitted to the hypernuclear proprieties. 
This can be done fixing the scalar a-meson hyperon coupling constants to the potential 
depth of the corresponding hyperon in saturated nuclear matter (U^ = —28 MeV, 
U£ = 30 MeV, US = -18 MeV) [51]. The obtained ratios x aY = g*Y 1 1 g<?N for the GM3 
parameter set are: x ct a = 0.606, = 0.328 and x aS = 0.322. Furthermore, following 
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Ref. [52], the SU(6) simple quark model can be used to obtain the relations 
1 1 1 

9 p n = ^9 P s = 9 P s , 9 p a = . (17) 

Concerning the formation of A-isobar matter at finite temperature and density, it 
has been predicted that a phase transition from nucleonic matter to A-excited nuclear 
matter can take place and the occurrence of this transition sensibly depends on the A- 
meson fields coupling constants [17] H8]. Whether metastable A-excited nuclear matter 
exists or not is still a controversial issue because little is actually known about the A- 
coupling constants with the scalar and vector mesons at finite temperature and density. 
On the other hand, QCD finite-density sum rule results predict a comparable or larger 
net attraction for a A-isobar than for a nucleon in the nuclear medium [19]. Because 
in this investigation we are going to focus on the nonextensive effects in hot and dense 
nuclear matter, we will not explore the existence of metastable A-matter and the scalar 
and vector coupling ratio will be fixed equal to nucleon one (g a A = g^N and g^A = g^N) 
(see, for example, Ref. [32] for more details about different values of A-coupling and 
the formation of A- metastable matter for different EOSs). 

The baryon pressure Pb = —Vt B /V and the energy density e can be obtained from 
the baryon grand potential Qb in the standard way. More explicitly, we have 

= | £ •» / S ww + "f(*)i - \ »4 *" - tf<-> 1 



2 2 

2 w 



+ ± c ( gujN uj) 4 + ^m 2 p p 2 , (18) 

f r\ 3 k 1 1 

^ = E * J 7^)3 K W + n?(k)] + -mla 2 + U(a) + - mW 

+ ^ c { 9ujN u) 4 + ^m 2 p p 2 . (19) 

The contribution of the lightest pseudoscalar (tt, K, K, i], t/) and vector mesons 
(p, u, K*, K* , 0) to the thermodynamics potential and to the others thermodynamics 
quantities, in the regime of temperatures and baryon densities explored in this work, 
is not negligible. On the other hand, the contribution of the ir mesons (and other 
pseudoscalar and pseudovector fields) vanishes at the mean- field level. Following Ref.s 
[3T| 132] . we can evaluate the pressure Pm, the energy density 6m and particle density 
Pm as that of a quasi-particle gas in the so-called one-body contribution, taking into 
account their effective chemical potentials depending on the self-consistent interaction 
between baryons. In this approximation we can write 
1 r r\ 3 k k 2 

p ^- 3 i:^Jw?wF) m+mi (20) 

= £ 7i / E M [jJW + , (21) 
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fi = n i^-Agm-Uk)]. (22) 



d 3 A; 

' 3 ,J J (2^p 

where jj = 2Jj + 1 is the degeneracy spin factor and the functions gj(k) and g~j(k) 
are the g-deformed particle and anti-particle distribution functions of the j-th meson. 
Again, for (3(Ej(k) — [£) > we have 

9j(k) = [l + (g-l)/3(^(fc)-^)]Vte-i) + i ' (23) 



where Ej{k) = w k 2 + m|. Moreover, the boson integrals are subjected to the constraint 
\f/,j\ < mj, otherwise Bose condensation becomes possible. 

The values of the effective meson chemical potentials fi* are obtained from 
the "bare" ones fij, given in Eq. (|T6|) . and subsequently expressed in terms of the 
corresponding effective baryon chemical potentials, respecting the strong interaction 
[3T1 132] . More explicitly, we have from Eq. ( TT6l) that /i 7r + = fi p + = \iq = p p — fi n and the 
corresponding effective chemical potential can be written as 

K+(p+) = Mp ~ AC — Pp ~ Pn - 9pN P , (24) 
where the last equivalence follows from Eq. (ll5j) . 

Analogously, by setting x^a = gw\/9uN, we have 

P*k+(k*+) = Pp ~~ Mace") = Pp ~ Ma — (1 — Xua)9ujnu — 2 9pNP , (25) 

Pk°{k*°) = Pn — A*A(s°) = Mn — Ma — (1 — Xu)a)9unV + 2 S'p^/ 5 > (^6) 

while the others strangeless neutral mesons have a vanishing chemical potential. Thus, 
the effective meson chemical potentials are coupled with the meson fields related to the 
interaction between baryons. This assumption represents a crucial feature in the EOS 
at finite density and temperature and can be seen somehow in analogy with the hadron 
resonance gas within the excluded-volume approximation. There the hadronic system 
is still regarded as an ideal gas but in the volume reduced by the volume occupied 
by constituents (usually assumed as a phenomenological model parameter), here we 
have a (quasifree) meson gas but with an effective chemical potential that contains the 
self-consistent interaction of the meson fields. 

Finally, the total pressure and energy density can be written as the sum of the 
baryon and meson contributions: P tot = Pb + Pm, ^tot = e_B + £m- 



3. Nonextensive QGP equation of state 

Concerning the nonextensive quark-gluon EOS, we use a simple effective MIT bag model 
to describe the quark-gluon phase in which all the non-perturbative effects are simulated 
by the bag constant B which represents the pressure of the vacuum. It is well known 
that, using the simplest version of the MIT bag model, at moderate temperatures the 
deconfinement transition takes place at very large densities if the bag pressure B is fixed 
to reproduce the critical temperature computed in lattice QCD. On the other hand there 
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are strong theoretical indications that at large densities (and not too large temperatures) 
diquark condensates can form, whose effect can be approximately taken into account by 
reducing the value of the effective bag constant [53] . A phenomenological approach can 
therefore be based on a bag constant depending on the baryon chemical potential [34] . 
In this study we consider a parametrization of the form 

B cS = (B - B 00 )/(l + exp[(/i fl - fio)/a]) + , (27) 

where we have fixed B^ 4 = 254 MeV (bag constant at vanishing fig), B^ 4 = 160 
MeV (bag constant at very large //#), /io = 600 MeV and a = 320 MeV. The above 
values have been obtained by requiring that, at vanishing chemical potential, the critical 
temperature is about 170 MeV for q — 1, as suggested by lattice calculations [54], while 
the other constraint is the requirement that the mixed phase starts forming at a density 
slightly exceeding 3 p f° r a temperature of the order of 100 MeV (as suggested by 
hydrodynamical calculations [55]). 

Following this line, the pressure, energy density and baryon number density for a 
relativistic Fermi gas of quarks can be written, respectively, as [50] 

d 3 A; k 2 
3 ^fj (2^~ f 

d 3 k 



_ poo r] 3 k k 2 

E/ o (28) 

_ poo j3l 

v E I (2^s e f + *?(*)] + B * . (29) 

(so) 

where the quark degeneracy factor for each flavor (/ = u,d,s) is 7/ = 6 and n,f(k), 
fif(k) are the g-deformed particle and antiparticle quark distributions. More explicitly, 

Uf{k) = [1 + (q ~ 1) P(e f (k) - »f)]W*-V + 1 ' (31) 



for (3(ef(k) — [if) > 0, where ey(fe) = J k 2 + m^. 

Light quarks m and d are considered as massless particles, while for strange quarks 
we consider a finite mass of m s = 150 MeV. Similar expressions for the pressure and the 
energy density can be written for gluons treating them as a massless g-deformed Bose 
gas with zero chemical potential. Therefore, the nonextensive pressure P g and energy 
density e g for gluons can be written as 

P^JlT — - (32) 

9 3 J (2tt) 3 [1 + (q - 1) fc/T]^-D - 1 ' 1 } 

e 9 =3P 9 , (33) 

with the gluon degeneracy factor j g = 16. In the limit q — > 1, one recovers the usual 
analytical expression: P g = (87r 2 /45) T 4 . 
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4. Mixed phase of hadron and quark-gluon matter 

As briefly mentioned in the Introduction, by increasing the baryon chemical potential, 
the phase diagram is characterized by a rapid crossover of thermodynamic observable 
around a quasi-critical temperature with a CEP, beyond which the system shows a 
first-order phase transition from confined to deconfined matter. 

The main goal of this study is to show the relevance of nonextensive statistical 
effects in a finite range of temperature and baryon density relevant for compressed 
baryonic matter experiments rather than the ultrarelativistic regime at vanishing baryon 
chemical potential. Thus, by considering the deconfinement transition at finite density, a 
mixed phase can be formed, which is typically described using the two separate equations 
of state: one for the hadronic phase and one for the quark-gluon phase. 

In order to study the phase transition from hadronic matter to QGP, we apply 
the Gibbs conditions to systems with more than one conserved charge, by requiring the 
global conservation of the baryon number, electric charge and zero net strangeness. The 
structure of the mixed phase is therefore obtained by imposing the following conditions 

(H) _ (Q) (H) _ (Q) (H) _ (Q) 

Hb ~ Hb ' He ~ He ' Hs — Hs ' 

PW(T, fJ, B , Hs) = P (Q \T, fj, B , fi C , fi s ) , 

PB = (l-X) Pb( T , Hb, He, Hs) + X Pb(T, [Ib, He, Hs) , ( 34 ) 

Pc = (1 - x) Pc ( T , Hb, He, Hs) + X Pc( T ' Hb, He, Hs) , 

Ps = (1 - X) Ps( t ,Hb,Hc,Hs) + xPs(T,Hb,Hc,Hs) , 

where p# , Pq^ and p§ are, respectively, the baryon, electric charge and strange 
densities in the hadronic (H) and in the quark (Q) phase and x is the volume fraction of 
quark-gluon matter in the mixed phase. At fixed T and the electric charge He and 
strangeness [is chemical potentials are obtained by fixing the net electric charge Z/A 
(for example, Z/A = 0.4 for lead-lead heavy ion collisions) and the net zero strangeness 
condition. 

As previously discussed, an important aspect it that the presence of more than one 
conserved charge implies a global and not local charge conservation, therefore the charge 
densities ps, Pc and ps are fixed only as long as the system remains in one of the two 
pure phases. In the mixed phase, the charge concentration in each of the regions of one 
phase or the other may be different. As we will see, this feature plays a crucial role on 
the strangeness production during the mixed phase of hadron and quark-gluon matter. 

5. Results and discussion 

Let us start our numerical investigation by reporting in figure [1] the variation of the 
pressure as a function of baryon density (in units of nuclear saturation density po = 0.153 
fm -3 ), at T = 120 MeV and for different values of q. Here and in the following we 
fix the value Z/A = 0.4. In presence of nonextensive statistical effects the pressure 
results to be considerably increased even for small deviations from the Boltzmann-Gibbs 
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Figure 1. Pressure as a function of baryon density (in units of the nuclear saturation 
density po) at T — 120 MeV and for different values of nonextensive parameter q. 

I [MeV] 




Figure 2. Phase diagram in the ps — T plane for different values of q. 

statistics. It is interesting to observe that the pressure presents a strong softening in the 
mixed phase even if the Gibbs conditions on the phase transition are applied. At this 
temperature, such a behavior results to be more pronounced by increasing the value of 
the nonextensive entropic parameter q with a larger range of baryon density involved 
in the mixed phase region. This matter of fact, already present in absence of strange 
particle degrees of freedom [29], implies an abrupt variation in the incompressibility 
and results much more evident here due to the additional strangeness conservation 
constraint. In this context let us observe that indirect indications of a remarkable 
softening of the EOS at energies reached at AGS have been already outlined [56] . 

In figure [2j we report the phase diagram in the ps — T plane for different values of 
the nonextensive parameter q. The curves labeled with p 1 ^ and p™ denote, respectively, 
the beginning and the end of the mixed phase. In presence of nonextensive statistical 
effects the phase diagram results significantly modified and a remarkable lowering of the 
critical maximum temperature at vanishing baryon density p_B is present. This result is 
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Figure 3. Variations of the baryon chemical potential pb (left panel) and the 
strangeness chemical potential ps (right panel) as a function of temperature at pu = Pa 
and for different values of q. 



in according to previous investigations where, by fitting the experimental observable at 
q > 1, the temperature (or slope) parameter T is usually less that the one obtained in 
the standard Boltzmann-Gibbs statistics (q = 1) j2ni EHJ 158] . 

In figure [3], we show the variation of the baryon chemical potential ps (left panel) 
and the strangeness chemical potential ps (right panel) as a function of temperature at 
fixed pB = Po and for different values of q. As expected, for a multicomposed strange 
hadronic matter, p s is positive and decreased with T at a fixed p B - It is interesting to 
observe the different behavior of ps at lower and higher temperatures in the presence 
(q 1) and in absence (q = 1) of nonextensive statistical effects. 

Let us now explore in more detail the particle concentrations for different values 
of the nonextensive parameter q. In figure IU we report the most relevant net particle 
ratios (Yj = pi/ Pb) as a function of baryon density at T=120 MeV. In presence of 
nonextensive statistical effects (right panel), the net particle concentrations are sensibly 
modified even for small deviations from the standard statistics. In particular, we observe 
in the hadronic phase a strong reduction of the neutron and proton fractions and a 
considerable increase in the hyperon and in the meson concentrations, also at moderate 
baryon densities. 

To better understand the relevance of the nonextensive statistical effects in presence 
of strange matter, we show in figure [5] the strangeness fraction Yg for baryons (B), 
mesons (M), strange quarks (s) and their antiparticles as a function of baryon density 
at a fixed temperature of T=120 MeV. For q = 1 (left panel), in the hadronic phase 
(x — 0), the total strangeness is carried almost completely by mesons (kaons, mainly 
K + and K°) and baryons (hyperons), although at low baryon density anti-mesons M 
(anti-kaons, mainly K~ and K ) bring a non-negligible fraction of strangeness. At 
the beginning of the mixed phase (p B ~ 3p ) the onset of s and s quarks rapidly 
dominates the strangeness ratio and the contribution of the other particle fractions 
becomes gradually less relevant. 

In presence of nonextensive statistical effects (right panel) the situation is quite 
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Figure 4. Particles concentration as a function of baryon density at T = 120 MeV 
and for q = 1 (left panel) and q = 1.1 (right panel). 



different. Here we observe a strong enhancement in all strange-particle concentrations 
and, in particular, the contribution of anti-baryons B (anti-hyperons) becomes relevant 
and comparable to anti-mesons and baryons. 




Figure 5. Strangeness fractions of baryons (B), mesons (M), strange quarks (s) and 
their antiparticles as a function of baryon density in the pure hadronic phase, mixed 
phase and quark phase at T=120 MeV for q — 1 (left panel) and q = 1.1 (right panel). 



This behavior is also evident in figure [6j where the strangeness fractions in the 
hadron and in the mixed phase as a function of the temperature at p = 3 po are shown 
(quark phase is not reported in this figure). As expected, by increasing the temperature 
all the strange particle densities are sensibly enhanced. Furthermore, in both panels we 
observe a peak at the end of the hadronic phase and a rapid decrease at the beginning 
of the mixed phase, due to the onset of the strange quarks. 

The strong enhancement of the strange hadronic particles in presence of 
nonextensive statistical effects is a direct consequence of the normalized mean 
occupation function because, for q > 1 and fixed baryon density (or /ig), it is 
enhanced at high values of its argument and depressed at low values. Being the 
argument of the mean occupation function Xi = j3(E* — p*), in the integration over 
momentum (energy), at lower ps (corresponding to lower values of p*) the enhanced 
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Figure 6. Strangeness fractions as a function of temperature in the pure hadronic 
and mixed phase at pb = 3 po f° r 9 = 1 (left panel) and q = 1.1 (right panel). 

nonextensive high energy tail weighs much more that at higher ps where depressed 
low energy effects prevail and the mean occupation number results to be bigger for the 
standard Fermi-Dirac statistics. Concerning the antiparticle contribution, the argument 
of Wi is Xi = (3(E* + p*) and the nonextensive enhancement at high energy tail is favored 
also at higher ps- As a consequence the formation of antiparticles results to be enhanced 
in presence of nonextensive statistical effects. At the same time, higher temperatures 
reduce the value of the argument of n« and Hi, reducing the effect of the nonextensive 
distribution function. 

These properties change significantly the behavior of the o meson field and, 
therefore, the effective baryon mass which is related by the relation: M* = Mi — g a iO~. 
In fact, for q > 1, M* becomes smaller at low baryon density (and higher temperature) 
and bigger at finite densities pb > 0.2 p (and lower temperature) [21]. This effect 
appears to be much more significant as a percentage for lighter baryons, therefore, at 
finite baryon densities, the nucleons effective mass is enhanced compared to the standard 
case (q = 1) with a percentage significantly greater than the hyperons effective masses. 
This matter of fact favors the formation of hyperons compared to that of nucleons at a 
fixed baryon density and in presence of nonextensive statistical effects. Moreover, being 
the strangeness number globally conserved, an increase of hyperon particles implies, in 
the pure hadron phase, a corresponding increase of strange meson particles in order to 
satisfy the condition of zero net strangeness. 

At the scope of better focalize the role of nonextensive statistical effects in the 
strangeness production during the mixed phase, in figure [3, we report the strange to 
anti-strange quark ratio p s /Ps a s a function of the baryon density, at different values 
of the volume fraction of quark matter \ in the quark-gluon phase along the phase 
transition, i.e., for a continuously varying temperature. In this context, let us point out 
again that, due to the Gibbs conditions, into the mixed phase, the net strangeness in 
each separate phase need not vanish although the total net strangeness is zero. This 
matter of fact is evident in figure [3 In the left panel (q = 1), especially at the beginning 
of the mixed phase (x = 0.1), there is a remarkable excess of quark s with respect to s, 
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with a maximum value around ps ~ 3 Pq. This means that in absence of nonextensive 
statistical effects, during the phase transition, there is a large antistrangeness (s) content 
in the hadron phase while quark-gluon plasma retains a large net strangeness (s) excess. 
This distillation mechanism, already known in literature, may result in "strangelet" 
formation, i.e. metastable droplets of strange-quark matter, which could imply a unique 
signature for quark-gluon plasma formation in relativistic heavy-ion collisions [57] . 

In presence of the nonextensive statistical effects (right panel), there is a very 
different behavior with a ratio p s j ' p-§ slightly less than one at lower baryon density 
(Pb < 2p ), whereas the other way round occurs at higher p B with a comparable 
behavior to the q = 1 case but with an excess of s quarks less pronounced. Similar 
behavior can be observed in figure where the same ratio as a function of the 
temperature is also reported. Of course, at fixed Xi l° w temperatures imply high 
densities and viceversa. 

As previous discussed, this behavior is a consequence of the power-law behavior 
of the mean occupational distribution for q > 1: the formation of antihyperons 
and kaons (mainly K + , K°) turns out to be disadvantaged at low baryon density 
and very high temperature and, viceversa, favored at high baryon density and 
low/intermediate temperature, compared to the standard case (q = 1). Although 
particles with strangeness and antistrangeness content are strongly enhanced in presence 
of nonextensive statistics, the ratio strangeness/antistrangeness results depleted due to 
an enhancement of antiparticles production at finite temperature. 

Ps/Ps ps/Ps 



1.5 




Figure 7. Ratio of strange to anti-strange quarks in the quark phase, as a function 
of the baryon density at different values of the volume fraction of quark-gluon matter 
X in the mixed phase for q = 1 (left panel) and q = 1.1 (right panel). 

In agreement with the previous results, in figure [9j we report the K + / K~ ratio in 
the hadronic phase as a function of baryon density at T = 120 MeV (left panel) and as 
a function of temperature at p B = 3p (right panel). For q = 1 and fixed temperature, 
the K + /K~ ratio increases with continuity by increasing the baryon density until it 
reaches the first critical density around (ps ~ 3 po), after that a strong enhancement of 
the ratio occurs in the mixed phase region. This a consequence of the antistrangeness 
(s) excess in the hadronic phase at the beginning of the mixed phase. For the same 
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Figure 9. Kaon to antikaon ratio as a function of baryon density (left panel) and 
temperature (right panel) for different values of q in the hadronic phase. 



reason the K + /K~ ratio decreases by increasing the temperature until the beginning 
of the mixed phase with an abrupt increase of the ratio. On the other hand, for q > 1, 
the kaon to anti-kaons ratio is depleted with respect to the standard case (q = 1), due 
to the increase of the antiparticles production. Furthermore, the strong variation of the 
ratio during the mixed phase transition at q = 1 results to be very smoothed in presence 
of nonextensive statistical effects. 

In presence of nonextensive statistical effects, the strong enhancement in density 
of strange particles compared to non-strange can be further observed in figure ITUl where 
the K + /tc + ratio is reported as a function of baryon density at T = 120 MeV (left panel) 
and as a function of temperature at p# = 3 po (right panel). For q = 1.1, there is a 
strong increase of the ratio, while it does not present discontinuities at the beginning 
of the mixed phase, unlike the case q — 1 where a much more antistrangeness excess is 
present. Finally it is interesting to observe for q = 1.1 a continuous decreasing trend of 
the ratio as a function of temperature at fixed baryon density. 
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Figure 10. The same of figure [9] for the K + /ir + ratio. 



6. Conclusions 

We have studied an effective nuclear EOS in the framework of nonextensive statistical 
effects at finite temperature and baryon dentity. By requiring the Gibbs conditions 
on the global conservation of the baryon number, electric charge fraction and zero 
net strangeness, we have studied the phase transition from hadronic matter to quark- 
gluon plasma and it has been shown that nonextensive statistical effects play a 
crucial role in the deconfinement phase transition also for small deviations from the 
standard Boltzmann-Gibbs statistics. Our investigation is focalized in regime of 
finite temperature and baryon density relevant for future compressed baryonic matter 
experiments. As previously discussed, in regime of high temperature and very small 
baryon density (low /ie), the phase transition may end in a second order critical endpoint 
with a smooth crossover and these features cannot be incorporated in the considered 
mean field approach. 

We have shown that strange particles are abundantly produced even at moderate 
temperature and the anti-particles concentration is strongly increased in presence of 
nonextensive statistical effects. For q > 1, the phase transition is characterized by 
an antistrangess content in the hadron phase while the QGP retains a net strangeness 
excess at large densities (p# > 2p ) and intermediate temperatures (T w 100 -v- 120 
MeV), while the other way round occurs at low densities and high temperature. This 
matter of fact is essentially due to the power law behavior of the mean occupational 
distribution function which weights differently low and high energy states at different 
baryon densities and temperatures. 

The possibility of separating strange from antistrange matter in the hadron-QGP 
phase transition can lead to a very significant enrichment of strange quarks in the QGP 
at high baryon density and intermediate temperature. In fact, in the hadronic sector 
of the mixed phase, the K + and K° are enhanced and the hyperons are suppressed. 
Being mesons much lighter than nucleons and their resonances, they carry away entropy, 
energy and antistrangeness and, therefore, the prompt kaon emission cools and charges 
the system with finite net strangeness, leading to an even stronger enhancement of the 
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s quarks in the quark phase. This feature, which also depends strongly on the value 
of the strange particle densities, favors the formation of metastable or stable droplets 
of strange quark matter which would contain approximately the same amount of u, 
d and s quarks. The evidence of such a state of quark matter could be related to 
the existence of exotic hadron states, like a if-dibaryon state, a deeply bound 6-quark 
state predicted by Jaffe more than thirty years ago [59]. Although many experimental 
searches for the if-dibaryon were carried out and so far no convincing signal was found 
[60] , very recently evidence for a bound if-dibaryon was claimed based on lattice QCD 
calculations [6T| 162]. 

In presence of nonextensive statistical effects, the separation of strange and 
antistrange quarks in the hadron-QGP mixed phase turns out to be less pronounced 
than in the standard case due to a more symmetric presence of particle and antiparticle 
at intermediate temperatures. On the other hand, for q > 1, the strangeness fractions 
Yg (and, consequently, the densities of strange particles) result to be much greater 
compared to the q = 1 case. This matter of fact could be crucial in the formation 
and survival of strange quark matter droplets in relativistic heavy ion collisions at high 
compressed baryonic matter. 
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